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We consider theoretically and numerically magnetic field dependencies of the maximum supercur- 
rent across Josephson tunnel junctions with spatially alternating critical current density. We find 
that two fiux-penetration fields and one-splinter-vortex equilibrium state exist in long junctions. 



I. INTRODUCTION 

Studies of periodic or almost periodic Josephson tun- 
nel structures arranged in sequences of interchanging 
0- and TT- biased Josephson junctions (as shown in 
Fig. [1]) recently became a subject of growing interest. 
These complex Josephson systems are intensively treated 
experimentally, theoretically, and numerically in: (a) 
superconductor-ferromagnet-superconductor (SFS) junc- 
tions in thin films^'^-'^''*'^ and (b) Josephson grain bound- 
aries in thin films of high-temperature cooper-oxide su- 
perconductor YBa2Cu307_:,,^ii^iaiiaiiiii^iHii4^ 

Equilibrium states of SFS Josephson junctions with a 
TT-shift in the phase difference between the superconduct- 
ing banks has been predicted almost three decades ago.^^ 
However, only recently SFS tt- shifted junctions and SFS 
heterostructures of interchanging - and tt - shifted frag- 
ments were studied experimentally for the first timci^'^'^ 

The studies of Josephson properties of the asymmet- 
ric grain boundaries in YBa2Cu307_2; thin films re- 
veal an interesting and important example of a Joseph- 
son system being an interchanging sequence of O-tt bi- 
ased j unctions i^ii'^i^'i^ The structure of these boundaries 
is created by facets with a variety of orientations and 
lengths / ~ 10— 100 nm.^- This grain boundary structure 
in conjunction with the (i^2_j,2-wave symmetry of the or- 
der parameter-^!^ can be considered as a Josephson tunnel 
junction with spatially alternating critical current den- 
sity, jc{x), where x axis is along the grain boundary^iSiiiS 
These rapid alternations with a typical length-scale of 
I significantly suppress the maximum supercurrent /,n 
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FIG. 1: Schematic drawings of: (a) sequence of 0- and tt- 
biased Josephson junctions; (b) critical current density, jc{x). 



across the grain boundaries. This suppression is most ef- 
fective for the asymmetric 45° [001]-tilt grain boundaries 
in YBazCuaOT-o. films-liiS 

The asymmetric 45° [001]-tilt grain boundaries in thin 
YBa2Cu307_a; films exhibit several remarkable and im- 
portant anomalies. First, the dependence of the maxi- 
mum supercurrent /„i on the applied magnetic field Ha is 
non-Fraunhofer.<2iiii^ii^ Contrary to the classical Fraun- 
hofer pattern with the central major peak two symmetric 
major side-peaks appear at the two fields ± Hsp ^ 0. Sec- 
ond, a spontaneous rapidly alternating magnetic flux is 
generated at the grain boundarieSf^ Third, unquantized 
spontaneous flux structures include fragments formed by 
pairs of single Josephson-type vortices carrying fluxes 
01 < 0o/2 and 02 > 0o/2'^^^ These fluxes are com- 
plimentary and sum to 0o, i.e., 4>i + 4'2 = 4>o and there- 
fore introduce splintered Josephson vortices. It is worth 
noting here that the anomalous patterns Im{Ha) and the 
unquantized splinter vortices appear under conditions of 
existence of equilibrium spontaneous flux. 

In many cases the length-scale I of the spatial alter- 
nations of the critical current density jc{x) is bigger or 
much bigger than the London penetration depth, A, and 
is smaller or much smaller than the local Josephson pen- 
etration depth, Aj, defined by the average of the abso- 
lute value of the critical current density. In the limit of 
/ ^ A J the phase difference between the banks of the tun- 
nel junction, (p(x), can be written as a sum of smooth, 
''pix), and rapidly varying, £,{x), termsJ^ Coarse-graining 
the phase ip{x) over a distance C » I allows to consider 
the two terms '((^{x) and ^(x) separately from each other 
in the inner part of the junction. In this approximation 
the coupling of 'ip{x) and S,{x) happens because of the 
boundary conditions at the edges of the junction. 

In this paper we calculate both theoretically and nu- 
merically the anomalous magnetic field dependence of the 
maximum supercurrent, /„, in Josephson tunnel junc- 
tions with spatially alternating critical current density. 
The applied magnetic field Ha is supposed to be lower 
than the side-peaks field, i.e., \Ha\ <S^ Hsp ^ 0o/27rAL 

The paper is organized as follows. In Sec. II we dis- 
cuss the coarse-grained equations for the phase difference 
across the banks of Josephson junctions with alternating 
critical current density and derive the boundary condi- 



tions to these equations. In Sec. Ill we consider the 
maximum supercurrent across Josephson junctions theo- 
retically in two limiting cases of short and long junctions 
in low and high magnetic fields. In Sec. IV we report 
on the results of numerical simulations of the maximum 
supercurrent dependence on the applied magnetic field. 
Sec. VII summarizes the overall conclusions. 



II. COARSE-GRAINED EQUATIONS 

We treat a one-dimensional Josephson junction par- 
allel to the X axis with the tunneling current density 
J II Yj jy{x) = j{x)i and the magnetic field H || z, Hz{x) — 
H{x). Assume also that the critical current density jc{x) 
is an alternating periodic or almost periodic function tak- 
ing positive and negative values with a typical length- 
scale I. The geometry of the problem is shown schemat- 
ically in Fig. [TJ 

First, we introduce the average value of the critical 
current density, ( jc), the effective Josephson penetration 
depth. A, defined by {jc), and the local Josephson pene- 
tration depth, Aj, defined by the average value of \jc\ 
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where Eq. ([T]) is the definition of averaging, L is the 
length of the junction (L ^ I), (Jjq is the flux quantum, 
and A is the London penetration depth. 

Next, we assume that A ^ Z <C Aj <C A. In this case 
the phase difference ip{x) satisfies the equation 
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It is convenient for the following analyses to write the 
critical current density jdx) in the form 



jc{x) ^ {jc) [I + g{x)] 



(5) 



introducing a rapidly alternating function g{x) with a 
zero average value, {g{x)) = 0, and a typical length-scale 
of order I. It is worth noting that g{x) is a unique internal 
characteristic of a junction. Using the function g{x) we 
rewrite Eq. dH) as 



AV- [l + g{x)] sinip^O. 



(6) 



The idea of the following calculation is based on a me- 
chanical analogy (Kapitza's pendulum) J^iii Two types 
of terms appear in Eq. ^ : fast terms alternating over a 
length I and smooth terms varying over a length A 3> /. 



The fast alternating terms cancel each other, indepen- 
dently of the smooth terms, which also cancel each other. 
Thus, to find solutions of Eq. ([6]) we use the ansatz 



ip{x) = ip{x) +^{x), 



(7) 



where ipix) is a smooth function with the length-scale 
of order A, ^(x) is a rapidly alternating function with 
the length-scale of order I, and the variations of ^{x) 
are small, i.e., {\£,ix)\) <S^ 1.^^ We assume also that the 
average value of ^{x) is zero, {^{x)) = 0. It is worth 
mentioning that the ansatz given by Eq. ([7]) is similar to 
the one used to solve the Kapitza's pendulum. ^^'^^ 

Substituting Eq. ([7]) into Eq. ([6]) and keeping terms up 
to first order in £^{x) we find^'^ 
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where the smooth j^{x) and alternating j^x) compo- 
nents of the tunneling current density j = j^ -\- j^ are 

Ji, = {jc) {sin ip- J sin il; cos -ip) , (10) 

Je = {Jc) g (x) sin ip, (11) 

the dimensionless constant 7 is equal to 

J = {g{x)^g{x)) , (12) 



and the rapidly alternating phase ^g{x) is defined by 



^{x) = -^g(a;)sin-0. 
It follows from Eqs. dH), ^ and ^ that 
AX+9ix)^0, 



(13) 



(14) 



i.e., the rapidly alternating phase shift ^g depends only 
on the effective penetration depth A and the function 
g(x). Therefore, the phase ^g{x) is an internal character- 
istics of a junction. 

It follows from Eq. (fTO)l that the smooth current den- 
sity j^ includes the initial first harmonic term ex simp 
and an additional second harmonic term oc sin 2ip, which 
results from constructive interference of the rapidly alter- 
nating critical current density ex g{x) and phase ^(x), — 

To summarize the derivation of the system of coarse- 
grained equations ((5|)- pT|) it is worth noting that the 
typical value of ^g{x) is small, but at the same time 
the typical value oi g{x) is big, i.e., {\^g{x)\) <C 1 and 
(|(7(a;)|) 3> 1. As a result, the dimensionless parameter 
7, which is proportional to the average of the product 
of the two rapidly alternating functions Cgix) and g{x) 
might be of the order of unity. ^■^•^^ 

The energy £ of a junction with alternating critical 
current density jdx) yields 



£ = 



2e 



„ (t*'^ 



^ ^" + 1 



cosTp — — sin^ V ) dx. (15) 



The last term in the integral in Eq. p5)) is for the con- 
tribution of both the fast alternating current j^{x) and 
phase S,{x). It is worth noting that minimization of the 
functional £{1^} results in Eq. ([5]) for the phase ipi^)- 

It follows from Eqs. ([8]) and (flS)) that if the parameter 
7 > 1, then there are two series of stable uniform equi- 
librium states with ipe = 2Trn ± ipj and current density 
jipi'ipe) = 0, where n = 0, ±1,±2, ... is an integer and 
the phase ^"7 is defined b}si^ 



7 cos ^^ — 1 . 



(16) 



All equilibrium states with ip — ip,, have the same energy 



_ h{j,) (7-1)^ , 



(17) 



which is less then the energy f = of the series of un- 
stable states with the phase ip = 2?™ J^ If the parameter 
7 < 1, then there is only one series of stable uniform 
equilibrium states with ipo = 2nn and £0 = 0. 

The two series of stable equilibrium states result in 
existence of two different single Josephson vortices (two 
splinters) M-^ The phase ip{x) for the first ("small") 



splinter vortex varies from —ipj at x 



-00 to ^py at 



X = -f cx). This vortex carries flux (pi = (pQip^/iT < (po/2. 
The phase for the second ("big") splinter vortex varies 
from ipj dA, X = ~cxd to 27r — -0^ at a; = cx3. This vor- 
tex caries flux (/)2 = 0o(7r — ■(/''y)/7r > 0o/2. As a result 
any flux structure inside a junction with an alternating 
critical current density and with 7 > 1 consists of series 
of interchanging small and big splinter vorticeSfi^i^^ It is 
also important mentioning that (pi + (p2 = 4'o- 

Consider now the boundary conditions to Eq. ([5]), i.e., 
for the smooth phase shift ip(x). Using equations 



H = _^rr, 



AttX dx 

ifi{x) = ipix) ~ ^g{x) smip{x) 



(18) 
(19) 



we find the boundary conditions for ip{x) in the form 

^'{0) ^ i;'o- ^so sin Vo = ^ ^0 , (20) 

ifi'iL) = ^P'l- ^gL sinV'L ^^Hl, (21) 

^ (po 

where Vo = ^(0), i^L = ^L), % = ^'(0), ^i = i^'{L), 
Ho = ff(0), Hl = H{L), C^o = e;(0), and ^^^ = Cg{L). 
Next, we use the fact that the average value of g{x) is 
zero and integrate Eq. (fT4| from to L. This leads to 



SeO ^ Cl — Seb' 
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where S.'^ is an internal parameter characterizing the 
edges of the junction. Now the boundary conditions given 
by Eqs. ^ and ^ take the form 



'fp'o - Cgb sin tpo 
^L-CgbSinV'L 
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Compare now the values of the derivatives ^'0, ip'i^, and 
Cgi,. Using Eqs. ^, ^ and dH]) we obtain 

and arrive to the relation 



A\ax)\^^^i. 



(26) 



A similar estimate A|?/''(a;)| ~ 1 follows from Eqs. ^ 
and pn]) . These estimates demonstrate that both deriva- 
tives ip'{x) and £.'g{x) are of the same order although 
(ICg(2;)|) ^ {\ip{x)\). Indeed, for a typical junction ex- 
hibiting spontaneous equilibrium flux we have 7 '^ 1.— 

The fact that A^'j^ ~ 1 makes it convenient for the 
following analysis to write the derivative ^' y^ in the form 



Li 



a 

A' 



(27) 



where a '--^ 1 is an internal parameter characterizing the 
edges of the junction. 

Thus, in the framework of the coarse-grained approach 
a junction with an alternating critical current density is 
characterized by two dimensionless parameters a and 7. 

Assume, that the current across a junction / 7^ 0, then 
we have the relations 



Ho ^ Ha^ / , 

c 



In this case the boundary conditions given by Eqs. 
and (12411 take the final form 



(28) 
(29) 
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(30) 
(31) 



The fact that the rapidly alternating critical current 
density jc{x) has low average value [{jc{x)) ^ (bc(2;)|)] 
might significantly affect the maximum supercurrent. In- 
deed, assume that the Josephson current density includes 
both the first and the second harmonicsF^ i.e.. 



j = jci (a;) sin (fi + jc2 sin 2(p , 



(32) 



where jci{x) is rapidly alternating along the junction and 
jc2 is spatially independent. 

In this case the coarse-graining approach remains the 
same as above. The effect of the second harmonics on the 
maximum supercurrent Im increases with the increase of 
the dimensionless parameter 72 = jc2/{jci)- The value 
of 72 might be of order of unity and higher even if jc2 is 
low compared to (|jci(a;)|)- 



III. MAXIMUM SUPERCURRENT 

The Josephson tunneling current, /, across Josephson 
tunnel junction with an alternating critical current den- 
sity can be written as a sum of two terms I^ and I^ 



I = jdx = I4, + /j , 
where the currents /^ and /^ are given by 



I^h — 



(33) 

(34) 
(35) 

(36) 



It follows from Eqs. (|M|) and ([55)1 that both /^ and /^ 
are defined by the smooth phase ip only. 
Magnetic flux inside the junction 







j^ dx^ IcA (iI^'l - -00 ) 



j^ dx ~ ale (sin -00 — sin-^i;,) , 



and the current I^ is defined as 

Ic = A{jc) 



= 2A / H dx = ^ i^L - t^o) 



results in the phase difference 



(37) 



(38) 



Using Eqs. ([55]) and ([55]) we obtain the current /j as a 
function of the flux inside the junction 



V'L - "00 = 27r -— . 

00 



le = ain 



sin -00 — sin ipQ + 27r^ 

00 



-2a/rSin ( tt — cos ■00 + ■"■ — 
00 / V 00 



(39) 



In order to calculate the total current I^ one has to 
know the spatial distribution of the phase ^^{x) in detail. 

In what follows we calculate the maximum super- 
current Im theoretically in the limiting cases of short 
(L <C A) and long (L ^ A) junctions treating the prob- 
lem separately for the Meissner and mixed states. 



A. Maximum current across short junctions 

We calculate now the maximum supercurrent I„i{Ha) 
of a short junction, L ^ A. In this case the spatial 
dependence of the smooth phase 0(a;) in the main ap- 
proximation in L/A ^ 1 is linear 



ip{x) = 00 + 27r 



(j>i X 



where 



2XLH, 



(40) 



(41) 



is the "internal" flux and Hi is the magnetic field inside 
the junction. Next, we use Eqs. dSO]), dST]), ^, (|35|) . 
and (|40p and obtain the following relations 



(42) 
(43) 
(44) 



Combining Eqs. (|43)) and (|44p we find that the maximum 
value of the total current /(0a) is given by 



Hi — Ha, 

0L = ^0 + 2^^, 4>a^2\LHa, 
00 

/ = ale (sin0L - sini/'o) • 



/„,, = 2alr, 



sm vr 



(45) 



It follows from Eq. (I45|) that the maximum supercur- 
rent across short junctions with spatially alternating crit- 
ical current density is defined only by the surface current 
/^ (in the main approximation in L/A ^1). As a re- 
sult the dependence Im{4'a) is obviously non-Fraunhofer. 
The value of /,„ is oscillating periodically in 0q with the 
period that is equal to flux quantum 00. Contrary to the 
case of a constant critical current density the amplitude 
of oscillations of Im is not decreasing with the increase 
of the apphed field HaM^ 



B. Meissner and mixed states in long junctions 

In this subsection we consider the spatial distributions 
of the phase difference and the flux in long junctions, 
L :> A. We start with the low field limit, i.e., we assume 
that the applied field Ha ^ Hs, where 



H,= 



2tt\A 



(46) 



is the flux penetration field for a long junction with a 
constant critical current density jc = const. ^^'^" In the 
following analysis we use an approach similar to the one, 
which was first developed by Owen and Scalapino.^i 

In the case of L ^ A and Ha ^ Hg the total super- 
current / = /^ + /j is a surface current localized in a 
layer with a width ~ Z ^ A. It follows from Eqs. ([34]) 
and (|35p that in order to calculate I^ and I^ we have to 
find the dependencies of V'i, and V'o on V'o and 0_l • These 
dependencies are given by the first integral of Eq. ([8]) 



A 



7 



— V -l-cosi/; cos 201 



const . 



(47) 



It is worth mentioning that Eq. ([T7|) describes the density 
of the energy £ given by Eq. p^ . 

The spatial distribution of 0(a;) depends on the mag- 
netic prehistory of the sample. We begin here for brevity 
with the case of a junction in the Meissner state. In this 
case the flux is localized at the edges of the junction. As 
a result in a long junction the phase ^-!{x) in the inner 
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FIG. 2: Spatial distributions of the phase "^(x) in a long junc- 
tion for different values of the applied field Ha and the internal 
flux <\>i. (a) Ha < Hsl, cj), = 0, (b) H,i < Ha < H,2, (t>^ = <^i, 
(c) Ha > Hs2 , (pi ~ po- 



part equals to a certain constant V'oo- The first correc- 
tion to this constant is proportional to exp (— L/A) <C 1. 
In other words we have 



1/^(^/2) = V'oo , i>'iL/2) = 0, 



(48) 



where the phase ipoo is given by one of the stable equi- 
librium values of ip, i.e., cosi/joo = 1/7- Combining Eqs. 
(|47)) . ((48)) and HH) we find that the constant in the RHS 
of Eq. (Ii7|) is given by 



const ~ cos t/iq 



(49) 



7 7 1 

-cos2V.^ = - + - 

The above relation allows for transforming Eq. (|47)l into 



A -0 = 7 (cos f/)^ — cos Tpy 



(50) 



We calculate first the flux penetration field into a junc- 
tion with an alternating critical current density and a 
zero total current, / = 0. The dependence ip^x) for this 
case is shown schematically in Fig.[5](a). It follows then 
from Eq. ^ that 



AV'o = V7(cos f/'o - cos ip^) , 
^''P'l — -\/7(cos'07 — cosiPl) ■ 



(51) 
(52) 

Next, we combine Eqs. dSO]), dST]), dSIl), and ^ and 
obtain two relations between the applied field Ha and 
the phases ■00 ^nd ■0^ 
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y/l + a^ cos (V'o + 0a) ^ 

V7. 



1 



VT + O^ cos (V^L - 0a) 



(53) 
(54) 



where we introduce the phase ipa as 

a 



tan il)a = — = 



(55) 



In the following analysis we assume, for definiteness, 
that "00- < ipj. In this case the dependence ip{x) looks as 
shown schematically in Fig. [51 

As a function of V'o the RHS of Eq. (JSH]) is bounded. 
The maximum field 



i/s 
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Vt^ 



1 



is achieved at 



■00 = -V'q • 



(56) 



(57) 



Therefore, if the applied field Ha reaches the value of 
Hsi then the Meissner state in a long junction becomes 
unstable and the small splinter vorte x^^i^"^ carrying flux 



io^ < 0o/2 



(58) 



enters into the inner part of the junction as shown in 
Fig. [21(b). This feature is a direct consequence of exis- 
tence of the splinter vortices in junctions with 7 > 1. 
It follows from Eq. (|50| that in this one-vortex state 



A0o,L = V7(cos0^ -cos-0o,l) > 0. (59) 

Using Eqs. pOj) . (|3T|) . and (|59|l we obtain the relations 
between the applied field Ha and the phases ipo and 0l 



Ha^^ 

2 



—- - \/'y + a^ cos (V'o - ipa) 
\/7 



Ha 



ff. 



V7 



VT + O^ cos (V-JL - 0a) 



(60) 
(61) 



The RHS's of Eqs. §0^ and ([61]) are bounded as func- 
tions of "00 and ipL said the maximum field 



ffs' 



2 



^7 + a^ -I 

\/7 



(62) 



is achieved at 0^0 = V'a — ti" + 27rn and ?/;£ = '>pa+TT + 2-Km 
where n,m = 0, ±1,±2,... are integers. If the applied 
field Ha reaches the value of Hs2 the one-vortex state 
becomes unstable and magnetic flux penetrates into the 
bulk until a mixed state with a finite density of vortices 
is cstabHshed (see Fig. [21(c)). 

Therefore, the rapid spatial alternations of the critical 
current density jc{x) in case of 7 > 1 lead to existence 
of a specific equilibrium one-splinter-vortex state. This 
state appear if the applied field Ha is from the interval 
Hsi < Hn < Hs2. It is worth noting here that the case of 
a standard Josephson junction (jc — const) corresponds 
to a = and 7=1. It follows then from Eqs. (|55|) . (|5^ . 
(|S^ . and (|^5)) that for these values of the parameters a 
and l3 we have 0i = 0, 02 = 00, ^si = 0, and i?s2 = Hs, 
i.e., there is only one Josephson vortex and the Meissner 
state exists if < Ha < Hg as it has to beJ^ This veri- 
fication means that the above results are self-consistent 
in describing the case of a standard Josephson junction. 



C. Maximum supercurrent in the Meissner state 

We calculate now the maximum supercurrent /„i in 
the Meissner state in a long junction, i.e., we assume 
that L ^ A and the smooth phase V' inside the junction 
is given by one of its equilibrium values ipe = 2?™ ± 
■07, where n = 0,±1,±2,... is an integer. The spatial 
distribution of ipi^) corresponding to the current /,„ is 
shown in Fig. [2](a). It follows then from Eq. (|5n|) that 



■00 = ^(cosV'^-cosV'o), 

iP'l ^ -jT- (COSV-L - COS ■07) . 



(63) 

(64) 



Using Eqs. ([55]) . ([M)) and the boundary conditions given 
by Eqs. ([50]) and pT|) we obtain equations relating the 
current /, the applied field Ha and the phases ipo and ^^ 



^a = 2 ^rn [cOs(V'L + ^a) - COs(-0O " V'a)] , (65) 



z 



Hr, 



2 

where we introduce the field Hm as 
iJ,i + Hs2 H, 

-Hm — z — — r 



COs(f/'o -1pa)+ COs(l/'L + l/^a)] , (66) 



\n + < 



(67) 



The two relations given by Eqs. (|65|) and (I66|) allow to 
obtain the dependence of the current / on the field Ha 
and the phase tpo in the form 



/c ^ Hs Hs 



(68) 



It follows from Eq. (|68p that the maximum current /„ 
corresponds to cos(-0o — ipa) = —1- Combining the above 
calculation valid for Ha > with the one valid for Ha < 
we obtain the dependence Im{Ha) in its final form 
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4/. 



H,- 
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27r 



iH,2-\Ha\)^ 



(69) 



Thus, in the Meissner state the maximum value of /„ 
is achieved at Ha = and is equal to 



/™(0) 



cH 



s2 



2tt 



= 21, 



V' 



7- 



1 ■ 



(70) 



It is worth noting that for a standard Josephson junc- 
tion (a = 0, 7 = 1) and therefore we have Hs2 = Hg. As a 
result Eqs. ([69|) and ([70| coincide with the similar equa- 
tions that were first derived by Owen and Scalapinoi^ 



D. Maximum supercurrent in the mixed state 

We calculate now the maximum supercurrent, /,„, in 
long junctions (L ^ A) in the mixed state, i.e., we as- 
sume that the applied magnetic field Ha is higher than 



i?s2- In the mixed state the field inside the junction, Hi, 
is almost uniform and tpix) takes the form 



^Hi X 

^ = V.o + 2 



(71) 



The dependence of the supercurrent on the applied field 
follows from the boundary conditions (|30p and (pij) yield- 
ing the system of equations 



4>a 4'i a L . f (pi 

'^~r = '^~l 77 T Sm-0m cos TT-— 

00 00 2 A V 00 



/ = — 2a/cCosV'mSin TT — 



where the phase </>,„ is defined as 

, -00 + V'L 

1pm = t; 



(72) 
(73) 

(74) 



Next, we use the Lagrange multipliers method to find 
the maximum of the supercurrent defined by Eq. (j73p 
under the constraint given by Eq. (j72p and arrive to 



ttIc 904 
00 dl 



C- 



Trie di)r, 



d(t)a 
di>m 



(75) 
(76) 



where C is the Lagrange multiplier to be determined. In 
the main approximation in A/L <^ 1 the solution of Eqs. 
((75)) and ([76]) is given by 



cosV'mCOS ( TT^^ I = ztsin^'mSin ( tt-^ 1 . (77) 

hj \ 00/ 



We plug now Eq. (177|) into Eq. (172|) and obtain 



_^ 27r A 0a - ( 
a L 00 



= cos 



(78) 



In the case of a long junction the LHS of Eq. (|78p is 
small. As a result in the zero approximation in A/L <C 1 
the flux inside the junction, 0^, is a constant defined by 
the roots of equation cos(7r0i/0o) = 0, i.e., the values of 
0i are given by 0^ = (n-l- l/2)0o, where n = 0,±1,±2, ... 
is an integer. In the next approximation in A/L <^ 1 the 
flux 0i depends on the flux 0^ and we find 



Ipn 




" + 2 ) 00 , (79) 



(80) 



where 



-ln+- 



(81) 



It follows therefore from the theoretical calculations 
that if the applied field Ha is increasing or decreasing, 
then inside the intervals [n — 1/2)00 < (j^a < {n + l/2)0o 
the flux in the bulk, 0,, is almost constant. At the ends 
of these intervals the flux (j)i "jumps" increasing or de- 
creasing its value by one flux quantum. 

Using Eq. (j73[) we find that the maximum supercurrent 
in the zero approximation in A/L ^ 1 is given by 

Im ~ 2a/, , (82) 

i.e., for long tunnel junctions {L 3> A) the value of /„ at 
high fields is almost field independent. 



IV. NUMERICAL SIMULATIONS 

We used numerical simulations to calculate the maxi- 
mum supercurrent in a wide range of parameters charac- 
terizing Josephson tunnel junctions with alternating crit- 
ical current density. The computations were performed 
by means of the time dependent sine-Gordon equation. 
The spatially alternating critical current density were in- 
troduced by the periodic function g{x). In the dimen- 
sionless form this equation yields 



Finite difference scheme 



(^ -H 5(^ - (^" + [1+5(0] sin </j = 0. 



(83) 



where the dimensionless time t — U,t and space C, — x/A 
variables are normalized by the Josephson frequency f2 
and length A, (5 <C 1 is the damping constant,— 



5(C) 



27rA 
27 — ; — sin 



'27rA 



V I 



■C + ' 



(84) 



the phase shift 9q defines the value of a, a = y/2^ cos Oq, 
and N = L/l is an integer {N ^1). 

The boundary conditions for Eq. ([83|) are given by the 
set of Eqs. ((281), ([291) and take the form 



(85) 
(86) 



^0 



Vl = 



2Ha 
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/ 
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2Ha 
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H, 




21. 



The convergency criterion for solutions matching equa- 
tions ([55]) and (|86p was based on the standard assumption 
that after sufficiently large interval of time (r ::^ 1) the 
spatial average of (^^(C,t) fits the condition {(p'^) < 5^, 
where (5m ^ 1 is a certain constant. We use a standard 
approach to calculate the maximum value of the super- 
current Im- Namely, for each value of the applied field 
Ha we find the current /,„ for which there is a solution 
of Eq. ([83ll matching boundary conditions ((85|) and ([86]) 
and converging after a certain time r, 3> 1, and there 
is no solutions converging at t ^ 1 for currents higher 
than /„j. We use the function ip{(,Tc) calculated for the 
field Ha as an initial condition <p((^, 0) for the next value 
of the field Ha + AHa, where AHa < Ha- 



We solved Eq. (|84p numerically using the leap frog 
method, which was adopted to our case. We checked the 
stability and convergency of the obtained solutions, and 
arrived at 
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(87) 
(88) 
(89) 
(90) 



where A^. and A^ are steps along r and C axes corre- 
spondingly, the superscript m is for the discrete r axis, 
and the subscript n is for the discrete ( axis. Next, we 
choose A^ to be equal to 1/12 of the period of the rapidly 
alternating function g{() and set A^ = A^. As a result 
we arrive at the following final difference scheme: 



V'' 



m+l 



= -il-6ArW:-' + i2-6A,)f>:: 
-A2(l+5„)sin(^™. 



(91) 



To obtain sufficiently accurate numerical data but to 
keep the time which is necessary for the numerical sim- 
ulations reasonable we choose the convergency criterion 
and the value of the decay constant S to be dependent 
on the length of the junction L. Specifically, we used for 
convergency criterion the following relations 



VT<^ < 10"^ for L < 8A, 
y/i^ < 10""*, for L > 8A . 



(92) 
(93) 



The value of the decay constant d of junctions with L < 
8A was chosen from S = 2 ior L = A/2 to S = .25 for 
L — 8A. In the case of junctions longer than 8A we took 
6 to be dependent on the convergency rate 



5 = 1.2, if ypy>io-^ 

6 = 0.1, if y/{^>io-^. 

B. Results of numerical calculations 



(94) 
(95) 



In this subsection we summarize the results of our nu- 
merical simulations for short {L <C A), long (L 3> A), 
and intermediate (L ^ A) junctions and compare the 
numerically calculated data to the theoretical results. 

In Fig. [3] (a) we demonstrate the dependence of the 
maximum supercurrent on the applied flux, Im{4'a), for 
a short junction, L = 0.25 A. In agreement with the the- 
oretical results obtained in Sec. IIII Al (see Eq. (|^ ) we 
find that Im{4'a) ~ Ic{4>a) except for small deviations 
at low fields. In Fig. [3] (b) we plot the internal fiux 




6 


*, ^ 


/ 




4 


/ 
/ 


J 
^ J 






r 


J 




" 


f 


J 







r __^ 


/ 
/ 


~^i 






J 


f 






J 


f 






J 


/ 






J ^- 


r 


a 


- 


^, ^\ 


/ 



FIG. 3: The maximum supercurrent Im and internal flux 
4>i dependencies on the apphed flux 0a for a short junction 
(L = 0.25 A, a = 2\/3 and 7 = 6). The arrows indicate the 
sweeping direction of the apphed flux, the points are for the 
results of the numerical calculations, and the solid lines are 
for the surface current lQ{4>a) given by Eq. (|45p . 
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FIG. 4: The maximum supercurrent Im and internal flux 0i 
dependencies on the applied flux </>„ for a short junction (L = 
0.5 A, a = 2\/3 and 7 = 6). (a) The function Im{4>a) with 
two side-picks located at (j>a ~ ±{L/l) c^o; (b) the dependence 
4>i[4>a) exhibiting flux-plateaus at the side-picks. 



4>i as a function of the applied flux (pa- It is seen from 
the graphs that (pi w 0a, which is in agreement with 
the assumptions of the theoretical calculations of Sec. 
nil Al Small flux "jumps", A(/) <C 0Oj are seen in Fig. [3] 
(b) in the vicinity of <j)a — ncjjQ, where n is an integer. 
These small flux jumps are generated by the high den- 
sity screening currents ^ (|jc(a;)|) ^ {jc{x)) flowing at 
the edges of the junctions. The length of these current- 
carrying edges is of the order of I and therefore the value 
of A0 can be estimated as follows. First, using Maxwell 
equations we find the field drop A7J at the edges to be 
Ai? « A-k{\Jc\)1/c. Next, we estimate At/) as a product of 
the field drop AiJ and the effective area of the junction, 
2AL, i.e., A0 « 2\LAH. Finally, we write the parameter 
q; as a « (|jc|) V(jc)A. Combining these three relations 
we find an estimate for A0 in the form 



A0« — -c 



(96) 



It is worth noting that A0 coincides with the coefficient 
in Eq. ([7^ for the difference between the internal flux 
and the applied flux. Using the data a = 2-^/3 and L = 
A/4 we obtain A0 = 0.1400 which is in a good agreement 
with the flux jumps shown in Fig. [3] (b). 

In this study we assume that the applied fleld is smaller 
than the side-picks field Hsp. The "resonances" at the 



side-picks are discussed in detail in Ref. '12 and 151 We 
show in Fig. [5] the maximum supercurrent Im{4'a) and 
internal flux (j)i{(j)a) at the side-picks for completeness and 
to reveal the flux-plateaus appearing in the dependence 

0i(0Q) at Ha = ±Hsp. 

In Fig. [5] we show the internal flux pi for a long junc- 
tion, L — BOA, as a function of the applied flux (jja- 
The value of (pi is less than one flux quanta if the field 
Ha is lower than the second penetration field Hs2- In 
this region of fields the slope dpi/ dp a is proportional to 
A/L ^ 1, i.e., it is almost zero. As a result, for long 
junctions in low applied fields we observe two relatively 
long flux-plateaus. These flux-plateaus, flux jumps and 
significant hysteresis in the magnetization curves pi{pa) 
are clearly seen in the whole area of pa- AH these fea- 
tures of magnetization curves are in a good agreement 
with the theoretical results obtained in Sec. IIIII 

In Figs. [S] (a), (b) we show the spatial distributions 
of the phase f{C) in a long junction, L = 30 A. The 
graph in Fig. [6] (a) is obtained for a junction in the 
Meissner state, i.e., for the applied field Ha from the in- 
terval < Ha < Hsi. In this case the flux inside the 
junction is zero. The graph shown in Fig. [5] (b) is cal- 
culated for a junction in the one-splinter-vortex interme- 
diate state, i.e., for the applied fleld from the interval 
Hsi < Ha < Hs2 and the internal flux pi = pi (see Eq. 
([55)) ). These numerical results are in a good agreement 
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FIG. 5: The internal flux (^,; dependence on the apphed flux 
<^a for a long junction (L = 30 A, / = 0, / = 0.2 A, a = and 
7 = 2). The arrows indicate the direction of sweeping of <^a 
for the magnetization cycles starting from t/ia = and zero 
trapped flux, (a) The main features of internal flux (jii{(j}a) 
of long junctions (flux-plateaus, flux jumps and significant 
hysteresis); (b) the dependence (j>i{4>a) for the applied field 
Ha from the interval < Ha < Hs2- 



with the theoretical calculation of Sec. IIIIBI 

In Fig. [7] (a) we plot the maximum value of the super- 
current Irn as a function of the applied flux (f>a for a long 
junction, L = 20 A. At low applied flux the maximum 
current is linearly dependent on \Ha\ yielding the middle 
triangle in agreement with Eq. (|69p . If the applied flux 
is sweeping up then the maximum current in the mixed 
state is higher than the maximum current in the Meissner 
state Ha = Hs2 — Analc/c and flux penetrates into the 
junction (flux jump) and the dependence Im{4>a) changes. 
If the applied field is sufficiently high then the maximum 
current is approximately equal to 2alc in agreement with 
Eq. ([H21). In Fig. [7] (b) we plot the internal flux <j)i as 
a function of the applied flux (pa- As it is assumed for 
fields lower than the first penetration field the junction 
is in the Meissner state. When sweeping the field from 
low to high fields the flux penetrates into the junction at 
Ha = Hs2 — inalc/c yielding a finite flux density. When 
sweeping the field from high to low values the Joseph- 
son vortices leave the junction one by one yielding the 
additional two steps between the plateau and the mixed 
state. In the interval of high applied fields the fiux jumps 
are of order of one flux quanta 4>q and in between the flux 
is almost constant in agreement with Eq. ([7^ . 

In Fig. [5] (a) we plot the maximum current as a func- 
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FIG. 6: The phase Lp dependence on the coordinate C, for a 
long junction (L = 30 A, J = 0, / = 0.2 A, a = 0, 7 = 2). 
The insets show the oscillatory nature of the function i-piC) on 
the space-scale of order /. (a) The applied field is sweeping 
up, the applied flux ^a = 0.8 0o; (b) the phase V'(C) is shown 
after the first fiux penetration that occurs at c^a = 4.6 4>o, the 
applied field is sweeping up. 
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FIG. 7: The maximum supercurrent Im and the internal fiux 
(jii dependencies on the applied fiux (j>a for a long junction 
(1/ = 20A, Q = \/3, 7 = 1.5). The arrows indicate the 
direction of sweeping of (j>a and the points are for the results 
of the numerical calculations. 
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FIG. 8: The maximum supercurrent Im and internal flux (^i 
dependencies on the applied flux (^^ for a junction with an in- 
termediate length (L = 2A, Q = 2%/3, 7 = 6). The points are 
for the results of numerical simulations. The arrows indicate 
the direction of sweeping of the applied flux. 



tion of the applied flux for a junction with an intermedi- 
ate length L = 2 A -- A. In Fig. [5] (b) we plot the flux 0,j 
as a function of the applied flux. It is seen that the flux 
0i differs from the applied flux 0a by less then one flux 
quantum (^q as for the short junctions. The flux jumps 



happen at (/)a = (n + 1/2)00, where n is an integer. The 
value of A0 is well approximated by Eq. (|M|) . 



V. SUMMARY 

To summarize, we consider theoretically and numeri- 
cally the maximum supercurrent across Josephson tunnel 
junctions with a critical current density, which is rapidly 
alternating along the junction. These complex Joseph- 
son tunnel systems were treated recently in asymmet- 
ric grain boundaries in thin films of high-temperature 
superconductor YBa2Cu307_j; and in superconductor- 
ferromagnet-superconductor heterostructures. 

Our theoretical study is based on coarse-grained sine- 
Gordon equation. We derive boundary conditions to this 
equation and find explicit dependencies of the maximum 
supercurrent across a junction on the magnetic field in 
the Meissner and mixed states for short and long junc- 
tions. We show that in the case of a Josephson junc- 
tion with rapidly alternating critical current density there 
can exist one-splinter-vortex mixed state and two flux- 
penetration flelds. The obtained theoretical results are 
verified by numerical simulations of exact sine-Gordon 
equation. We demonstrate that the theoretical and nu- 
merical results are in a good agreement. 
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